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The quantum-classical transition of the escape rate of the spin model TL = —DS Z — H Z S Z + BS X 
is investigated by a perturbative approach with respect to B [D. A. Garanin, J. Phys A 24, L61 
(1991)]. The transition is first order for B < B C (H Z ), the boundary line going to zero as B c /D ~ 
1 — H Z /(2SD) in the strongly biased limit. The range of the first-order transition is thus larger 
than for the model TL = —DS Z — H Z S Z — H X S X studied earlier, where in the strongly biased case 
H XC /(2SD) ~ [1 — H z /(2SD)] 3l/2 . The temperature of the quantum-classical transition, To, behaves 
linearly in the strongly biased case for both models: To ~ 2SD — H z . 

PACS number(s): 75.45. +j, 75.50.Tt 



I. INTRODUCTION 

Recently we have shownfli that the order of the 
quantum-classical escape-rate transitionu in the spin sys- 
tem described by the Hamiltonian TL = —DS Z — H X S X 
is controlled by the transverse field H x . According to 
Ref. H for the reduced field h x = H X /(2SD) in the 
interval 1/4 < h x < 1 the transition is second or- 
der, which is a common situation (for h x > 1 the bar- 
rier between the two wells disappears). In contrast, for 
< h x < 1/4 there is a first-order transition of the 
escape rate T(T) characterized by the discontinuity of 
dY/dT at the transition temperature To in the large- 
spin limit. Subsequent calculations^ rendered the whole 
phase diagram for the model with an arbitrarily directed 
field, TL = -DS 2 - H Z S Z - H X S X , where the first-order 
transition occurs for h x below the line h xc (h z ), where 
h z = H Z /(2SD). This line starts from the value 1/4 
at h z = and approaches zero as h xc ~ (1 — h z ) 3 / 2 
in the strongly biased limit. This limit can be rele- 
vant for observation of the crossover between first- and 
second-order transitions on single-domain magnetic par- 
ticles, where the barrier should be lowered by applying 
the field to make escape rates measurable. On the other 
hand, for molecular magnets, such as Mni2 (S = 10, 
D = 0.6 K), experiments can be done even in the unbi- 
ased case, H z = 0. 

The mechanism leading to an exotic first-order escape- 
rate transition in spin systems is the following. In the ab- 
sense of the transverse field the eigenstates of the system 
are those of the operator S z , and there is no tunneling 
between the wells. This implies that the barrier of the 
effective potential of the spin system, which can be ob- 



tained by a mapping onto a particle problemJS'B becomes 
infinitely thick in the limit H x — > 0,preserving, however, 
its height AU (see Fig. 1 of Ref. |). That is, the bar- 
rier has a very flat top and it resembles a rectangular 
barrier. Tunneling just below the top of such a barrier 
is very unprobable. The thermally assisted tunneling is 
thus suppressed, and the thermal activation competes di- 
rectly with the ground-state tunneling, leading to a sharp 
escape-rate transition. 

Mapping of a spin problem onto a particle one is, how- 
ever, not a regular procedure, and its form strongly de- 
pends on the form of the spin Hamiltonian (see Ref. ||). 
In many cases, such as for the model with S z , it is dif- 
ficult to find a particle mapping. On the other hand, 
one needs a simple criterion for roughly estimating the 
order of the transition for different spin systems. This 
criterion can be constructed as follows. The first-order 
transitions occur if a spin system is close to a uniaxial 
one, TL = F(S Z ). In this case the maximum of the clas- 
sical energy of the spin, U(S), corresponds to a certain 
value of the polar angle, 9 — 6 C , and there is no saddle 
point of U(8, ip). If nonuniaxial terms, such as the trans- 
verse field, are added to the Hamiltonian, then a saddle 
point appears. This brings the spin system closer to the 
common case, and if the saddle is strongly pronounced, 
one can expect a second-order transition. Thus, one can 
propose the following heuristic criterion for the boundary 
between first- and second-order transitions: the depth of 
the saddle is of the order of the height of the barrier, i.e., 
t^max — £4ad ~ C/ Sa d — U nlln . For the model with the trans- 
verse field one has AU = U sa d — U m i n = S 2 D(1 — h x ) 2 
and t/max — U s& & = AS 2 Dh x . Equating these expressions 
one obtains h xc = 3 — 2\/2 ~ 0.17, which has a proper or- 
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der of magnitude. In the case H z ^ 0, the parameters of 
the barrier cannot be calculated analytically, but estima- 
tions can be done. In particular, in the strongly-biased 
case, 6 = 1 — h z -C 1, for h x = the top of the barrier 
corresponds to 9 C = n — 9 C = V%8, and AU = S 2 DS 2 . 
For h x ^ 0, one can use U max — U sa d ~ S 2 Dh x 9 c and the 
same barrier height for estimations, which results in the 
correct dependence h xc ~ <5 3 / 2 (Ref. |]). 

In this communication we will study the quantum- 
classical escape-rate transition for the spin model with 
a biaxial anisotropy, which is described by the Hamilto- 
nian 

H = -DS 2 Z -H Z S Z + BS 2 X . (1) 

This model has the easy axis z, hard axis x and thus 
the middle axis y. Applying the longitudinal field is the 
simplest way to reduce the barrier. Sometimes another 

form of the biaxial spin model, TC = K(S 2 + XS 2 ) H , 

is used, which for A < 1 has the easy axis x, hard axis 
z and the middle axis y. These two forms are related by 
D = \K and B = (1 — \)K. A realization of this spin 
model is the magnetic molecule Fes (S = 10, D = 0.31 K, 
B = 0.092 K, Refs. § and @). 

In the unbiased case, H z = 0, using the criterion 
C^max - C4ad = U s&d - C^min readily yields b c = 1 for 
the critical value of the reduced quantity b = B / D corre- 
sponding to the boundary between the first- and second- 
order transitions. Occasionally, this coincides with the 
exact value of b c which has been obtained recently in Ref. 
|8|by the periodic instanton method (A c = 1/2). In the bi- 
ased case, for the reduced energy U = U /(S 2 D) one finds 
Umin = -1 + 2h z , U sad = h 2 z , and L W = b + h 2 z /(l + b). 
Our criterion yields b c = (1 — h z )(yjl + h 2 — h z ), i.e., b c 
decreases linearly with 8 = 1 — h z in the strongly biased 
case. Thus, in this case the region where a first-order 
transition can be expected, is wider for the biaxial model 
than for the model with a transverse field (h xc ~ <5 3 / 2 ). 

Now we proceed to an actual calculation of the bound- 
ary line b c (h z ) for the model above. In principle, this 
could be done by a kind of instanton or other quasiclas- 
sical method. However, the problem of obtaining rele- 
vant instantons (or even the action without computing 
instantons explicitly) in the whole field range appears 
to be mathematically difficult. Thus we choose to solve 
the problem with the help of the high-order perturbation 
theory with respect to b = B/D in Eq. ([!]). The lat- 
ter was applied to calculate the splittings of the ground 
and excited states for the unbiased model in Ref. ||. It 
was shown that the instanton results of Ref. [l(] for the 
ground-state splitting are recovered in the limit ii « 1, 
and that the Kramers degeneracy for half-integer values 
of S appears in a natural way. For the model with a 
transverse field, this method was used in Ref. [ll] to cal- 
culate ground-state splitting in rare-earth compounds. It 
has been generalized for excited states in Ref. ^. Similar 
results have been obtained in Ref. |l2| for the model of two 
coupled nonlinear oscillators (quantum dimer), which can 
be considered as a Schwinger-boson equivalent of the spin 
model with a transverse field. 



The perturbative method was first to show the first- 
and second-order transitions in the spin model with a 
transverse field in Ref. [jj. The perturbatively determined 
value h xc — 0.145 was, however, wrong by a numeri- 
cal factor in comparison to the exact value h xc ij 1/4 
obtained subsequently by a quasiclassical method.EI The 
reason for this inaccuracy is that the order of the escape- 
rate transition is controlled by the situation near the top 
of the barrier, just where the perturbation theory breaks 
down. A plausible way to improve the accuracy of the 
perturbative results is to rescale them by fitting to the 
exact quasiclassical results at h z = 0. Then one can ex- 
pect improving the results in the whole range < h z < 1, 
where, for the biaxial model, exact quasiclassical calcu- 
lations have not yet been performed. At first we will 
test this method on the model with a transverse field, for 
which the quasiclassical boundary line h xc (h z ) is known 
from Ref. [|. We will see that, surprisingly, a simple 
rescaling of the perturbative results yields the accurate 
quasiclassical ones in the whole range < h z < 1. 

Within the perturbative method, the quantum- 
classical transitions correspond to different types of be- 
havior of the functiontl 

/(m) = (Ae mm ,) 2 exp(-e m /T), (2) 

where £ rn = —Dm 2 — H z m are unperturbed energy lev- 
els of the spin system and Ae mm > is the splitting of the 
pair of in-resonance levels e m and e m < belonging to dif- 
ferent potential wells. This resonance condition is ful- 
filled for the values of the bias field H z = H z k = kD, 
where fc = 0, ±1,±2, . . .; then one has vn! = — m — k. 
The function /(m) is the quantum- mechanical tunneling 
probability T mm i cx (Ae mm ') 2 weighed with the Boltz- 
mann exponent exp(— e m /T) characterizing the proba- 
bility of thermal activation to the level m. The task is 
to find the escape level m OS c, i-e., the value of m which 
maximizes /. The maximum of / is searched at the in- 
terval — S<m< TTtjj , where m — —S corresponds to 
the bottom of the well and m = mf, corresponds to the 
top of the barrier which is determined from the condi- 
tion that the level splitting Ae mm ' reaches the value of 
the level separation e m — e m _i (see, e.g., Ref. Q). For 
temperatures above the transition temperature To, the 
result will be m cac = mj,, which corresponds to thermal 
activation over the barrier. For T < To, one obtains 
m esc (T) < nib, which means thermally assisted tunnel- 
ing (or ground-state tunneling for m csc = — S). If the 
dependence m csc (T) is smooth, the spin system under- 
goes a second-order escape-rate transition. If there is a 
jump of m esc (T) at some temperature, there is a first- 
order transition. 

For the model with an arbitrarily directed field and 
without transverce anisotropy the level splittings are 
given by the formula 

Ac , — 

[{m — m — 1 

l (S + m>y.(S-m)\ (H x \ m '- m 
V {S-m')\{S + m)\ \2D) ' [ ) 
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FIG. 1. The level making the dominant contribution into 
the escape rate m eS c vs temperature for the unbiased trans- 
versed field model. 
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FIG. 3. Quantum-classical transition temperature To at 
the boundary between first- and second-order transitions for 
the model H = -DS 2 Z - H Z S Z - H X S X . 
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FIG. 2. Boundary between the first- and the second-order 
transitions for the model TC — —DS 2 Z — H Z S Z — H X S X . 
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FIG. 4. Boundary between the first- and the second-order 
transitions for the model Ti. = —DS Z — H Z S Z + BS X . 



which is a generalization of the zero-bias result of Ref. g 
(see, e.g., Ref. |l|). Numerical maximization of the func- 
tion /(to) of Eq. (||) leads to the results for m csc (T) 
shown in Fig. [l]. One can see that the transition is sec- 
ond order for h x > 0.145 and first order for h x < 0.145. 
The results in Fig. [l] suggest that in some range of h x 
(as, e.g., for h x = 0.13) there are two transition temper- 
atures: With lowering temperature, TO esc at first devi- 
ates from TOb (a second-order transition) and then jumps 
downwards (a first-order transition). In fact, this feature 
is an artifact of the perturbative approach to tunneling, 
which breaks down near the top of the barrier. (For this 
reason also the crossover field h xc = 0.145 substantially 
deviates from the exact value h xc = 0.25 obtained in Ref. 
||.) In the following we will associate the crossover from 
first- to second-order transitions with the transverse field 
at which the jump in the dependence m esc (T) appears. 

The values of the crossover field h xc for different lon- 
gitudinal fields h z , which have been obtained by the 
perturbative method described above, are shown by the 



solid symbols in Fig. |2j. These data can be corrected 
by multiplying by the factor 0.25/0.145=1.274 to fit to 
the exact quasiclassical value h xc — 0.25 in the unbi- 
ased case. Surprizingly, this rescaling leads to the results 
which completely coincide in the whole range of the bias 
field with the results obtained in Ref. [| by the quasiclas- 
sical method based on the particle mapping. 

The results for the transition temperature To at the 
boundary between first- and second-order escape-rate 
transitions in the whole range of the bias field is shown in 
Fig. ^. Again, the perturbative results can be corrected 
by multiplying by 0.1378/0.1198=1.150 to fit to the exact 
value in the unbiased case: To(h xc ) = [\/3 / (Att)]S D = 
0.1378 SD (Ref. ||). This makes them accurate in the 
whole range of h z , as shows the comparison with the ex- 
act asymptote T (h xc ) ^ 0.1642(1 - h z )SD (Ref. |) in 
the strongly biased limit. 

Having tested the perturbative method with rescaling 
on the model with an arbitrarily directed field, let us now 
proceed to the biaxial model of Eq. [fl, where a quasiclas- 
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FIG. 5. Quantum-classical transition temperature To at 
the boundary between first- and second-order transitions for 
the model H = -DS 2 Z - H Z S Z + BS 2 X . 



sical solution is still lacking in the biased case, H z ^ 0. 
For this model the formula for the level splittings reads 



As, 



8D 



{[(m'-m)/2- l]!} 2 

' (S + m'y.jS-m)] f B \ ( m '-"*)/ 2 
{S - m')\(S + m)\ \\QDJ 



(4) 



less justified, and one can expect some deviations of the 
rescaled perturbative results from the exact quasiclassi- 
cal ones for H z ^ 0. Thus obtaining of the latter is an 
actual problem. 

For single-domain magnetic particles with small trans- 
verse anisotropy B, i.e., for nearly uniaxial ones, the 
crossover from first- to second-order escape-rate tran- 
sitions should be searched for in the strongly biased 
case, 5 = 1 — h z <C 1, where the barrier is reduced. 
Here our result for the transition temperature at the 
boundary between first- and second-order transitions is 
T (B C ) « 0.39 SDS (see Fig. |). On the other hand, 
the unperturbed barrier height is given by AU — S 2 DS 2 . 
Eliminating S from the above formulas, one can relate 
T and AU as follows: T « (0.39) 2 T>A£//T o . The ratio 
AU /To is fixed by the requirement that the escape rate 
be not too high and not too low. Adopting a typical value 
AU /To = 40, as was done in Ref. ^ for the model with 
an arbitrarily directed field, one arrives at the estimation 
Tq(B c ) « 6. 3-D, which is substantially higher than that 
of Ref. |, T (h xc ) D. 
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ogy with the quantum dimer problem and for the critical 
reading of the manuscript. 



which is a generalization of the zero-bias result of Ref. 
^. This formula is only applicable if m' — m is even, 
otherwise the perturbative method yields As mm i = 0, 
i.e., the Kramers degeneracy. Hence, for S integer and 
the resonant values of the bias field H z = H z k = kD, 
the level pairs are splitted for k — 0, ±2, ±4, . . ., but re- 
main degenerate for k = ±1,±3, . . .. For S half integer, 
the situation is reversed. Below we will only consider 
the resonant values of H z for which there is no Kramers 
degeneracy. 

The dependences m esc (T) for the biaxial model are 
similar to those for the model with an arbitrarily di- 
rected field, and they are assessed in the same way. The 
results for the bondary between first- and second-order 
transitions, B c , are shown in Fig. ||. For rescaling of 
these results we have used the exact value B c /D = 1 
at H z = (Ref. ||). One can see that B c cx 1 — h z 
in the strongly biased case, as was conjectured above. 
The transition temperature To at the boundary between 
first- and second-order transitions in the whole range 
of the bias field is shown in Fig. ||. Here we have 
used for the rescaling of the results the exact value 
T (B C ) = (l/n)SD = 0.3183 SD at H z = 0, which can 
be extracted from Ref. ||. 

In comparison to the previously considered model with 
an arbitrarily directed field, here the crossover between 
first- and second-order transitions occurs at higher val- 
ues of the control parameter B, thus the range of the 
first-order transitions is extended. On the other hand, 
for this reason using a perturbation theory in B/D is 
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